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Abstract 

A stochastic heat equation on an unbounded nested fractal driven by a general stochastic 
measure is investigated. Existence, uniqueness and continuity of the mild solution are proved 
provided that the spectral dimension of the fractal is less than 4/3. 

1 Introduction 

The main object of the paper is the stochastic heat equation on an unbounded nested fractal E 
which can formally be written as 



measure, C is the infinitesimal generator of Brownian motion on E. We consider the equation in 
the mild sense (see ([7]) below) and obtain existence, uniqueness and continuity of the solution. 
The Euclidean version on the real line was proved in [10J. 

Here we use the Brownian motion on E constructed in Our results are based on the Holder 
continuity and an upper estimate of the transition density of the process derived in that paper. 
Diffusions on nested fractals also were constructed in [TJ EJ [6]. 

SPDE on fractals and general metric measure spaces are studied e.g. in |7j (with an infinite- 
dimensional Wiener process as integrator) and in [2] (with pathwise definition of the stochastic 
integral). In our paper we consider the case of an additive noise, but with a rather general stochastic 
integrator. We do not assume path regularity or moment existence for the stochastic measures. 

The paper is organized as follows. Section [2] summarizes some basic facts about nested fractals 
and stochastic measures. In Section [3] we obtain the continuity of paths of stochastic parameter 
integrals on these fractals. Section H] contains the main result of the paper. 
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2 Preliminaries 



2.1 Heat kernel on nested fractals 

Let E be a nested fractal (see Definition 1.3 [4j). Then E C IR d is a compact set, E has a self-similar 
structure, i.e., E = {J^ip^E), where ^ : R ' — >■ IR d are similitudes with the same contraction 
factor, 

\ipi{x) - ipi(y)\ = a~ x \x -y\, a> 1. 

In the following we assume that ip\(x) = a~ l x. 

Let us denote by d s the spectral dimension, by d w the walk dimension, and by df the Hausdorff 
dimension of the fractal set E. Recall that df = log N/ log a, = 2df/d s . 

Let F be the set of fixed points of 1 < i < N. x G F is called an essential fixed point if 
there exist y G F, j ^ k such that ipj( x ) = ^k{y)- Let F^ be the set of essential fixed points, 

F^= |J ^o..-o^(F(°)). 

Further, let us define the unbounded nested fractal E. We set 

= |J (a^W), = S = CI ((J fW) . 

n>0 nSZ 

Note that 

£=U U K^o-.-o^^)), (2) 

n>l ii,...,i n 

where any two sets in the union either coincide or have finite intersection. 

Let m be the df- dimensional Hausdorff measure on E such that m(E) = 1, 13(E) denote the 
Borel cr-algebra on E. 

Assumption 1. There exists fceN satisfying the following. If for x, y G E we have \x — y\ < 
a~ m , then there exist xi, X2, ■ ■ ■ , xi (I < k) such that x\ — x, x\ = y, x%, . . . , G F^ m ^ and 
Xj, 1 < j < I — 1 lie in the same set ip^ o • • • o tej . 

Let p(t,x,y), t > 0, x, y G 22, be the transition density of Brownian motion on E constructed 
in |4J provided that Assumption 1 holds. Then p is symmetric in x, y, is jointly continuous in 
(t, x, y) and satisfies the following Holder continuity and sub-Gaussian estimate for some constants 
ci, C2, C3 > 0, and dj > 1. 

Theorem 1. Theorems 5.1 and 5.2] 

a) \p(t,x,yi) -p(t,x,y 2 )\ < cxt^yx -y 2 \ / . 

b) p(t,x,y) < c 2 r^/ 2 exp{-c 3 (|x-y| d -/t) 1/(dj_1) }. 

We will take this p(t, x, y) as the heat kernel for the stochastic heat equation in the mild 
form ([7]). Note that p fulfills the deterministic heat equation dp/dt = Cp, where £ is the infinites- 
imal generator of the Brownian motion on E [U Theorem 6.3]. 

The same estimates of the heat kernel on the Sierpihski gasket, both the bounded and un- 
bounded cases, are given in Theorem 2.23 [lj, and the inequalities for the two cases are equivalent. 
In our paper, we consider unbounded nested fractals, but we expect the same results for the 
bounded versions. 
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Example 1. The Vicsek set is the main example of this paper. Let H be the closed unit square 
with vertices a\ = (0, 0), a2 = (0, 1), a% = (1, 1), = (1, 0). We take = (1/2, 1/2) and set 

5 

ipi{x) = di + (x - Oi)/3, 1 < i < 5, H n+ i = (J ipi(H n ), E vs = f] H n . 

i=l n>0 

Then Eys is called a Vicsek set (or Vicsek snowflake). 

It is known that Eys is a nested fractal |3j Proposition 2.1], and we have for this set df = 
log5/log3, d w = logl5/log3 [TJ Section 2], d s = log25/logl5 (see calculations in (TTJ Section 6] 
for n = 2). Note that Assumption 1 holds for this set with k = 4. 

2.2 Stochastic measures 

Let L = L (Q, J-, P) be the set of all real- valued random variables defined on the probability 
space (Q, J 7 , P) (more precisely, the set of equivalence classes). Convergence in Lq means the 
convergence in probability. Let X be an arbitrary set and B be a u-algebra of subsets of X. 

Definition. A a -additive mapping fi : B — > L is called stochastic measure (SM). 

We do not assume any martingale or moment conditions for fi. Examples of SM are the fol- 
lowing. For square integrable martingale X(t), < t < T, set fi(A) = Jj T , Ia(^) dX(t), then 

fi is a SM on the Borel cx-algebra B([0, T]). For unconditionally convergent in probability series 
E„>i an d {x n } C X set function fi(A) = X) n >i ^A{x n ) is a SM on 2 X . 

Example 2. Let Eys be the Vicsek set. We construct an example of a SM on the Borel <r-algebra 
B[Eys\ Let /i be a SM on B((0, 1]) such that /i({x}) = a. s. for each x G (0, 1]. For n > we 
take a representation 

Evs = u 4" ) =^o---o^(^)=^i ) u---u^r i) . 

l<fc<5 n 

We set 

IMe(eM) = ^((k - l)5-\ A;5-]), ^({x}) = 

for each x G C\E\ and define /x^ on the generated algebra by additivity. Then fig is continuous 
in and the values of \ie are bounded in probability (see (5J Theorem B.2.1]). By Theorem 1 of j8] 
we can extend the stochastic set function to B{Eys)- Taking the weighted sum of stochastic 
measures defined on the sets a^ip^ o • ■ • o ip in [Eys) from (j2J), we can easily extend this stochastic 
measure to the unbounded Vicsek set Eys- 

In [5J the integral of measurable deterministic functions / : X — » R with respect to (w.r.t.) SM 
is constructed. The dominated convergence theorem holds for this integral [5J Proposition 7.1.1], 
and all bounded measurable functions are integrable. We will use the following statement. 

Lemma 1. [10, Lemma 3.1] Let g\ : X — >■ R, I > 1, be measurable functions such that g(x) = 
i s integrable w.r.t. jj,. Then Yl'iZiilxd 1 ^^) < 00 a - s - 
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3 Parameter stochastic integrals on nested fractal 

Let E be an unbounded nested fractal, /x be a SM defined on 13(E), Z be a metric space, function 
h : Z x E — > K be integrable w.r.t. /ion B for each first argument 2 G Z. Then the stochastic 
function 

77(2) = / dfi(y), z G Z 

is determined. 

Theorem 2. Suppose that h is bounded and continuous on Z x E, for some Kh > 0, (3(h) > df/2 
we have 

\h(z,y x ) - h(z,y 2 )\ < K h \y x - y 2 \ m , (3) 

and the SM /j is such that fi({x}) = a. s. for each x G E, and for some r > df/2 the function 

\y\ T is integrable w.r.t. /j on E. 

Then rj(z) has a version with continuous paths. 

Proof. Using (J2J), we can take a representation E = Uj>i Ej, where each Ej is a bounded part of 
the form 

a> n o ■ ■ ■ o t/j in (E), Ej C a n E, Ej <£ a n ~ x E for N 71 ' 1 <j< N n . 

Let 

7 = inf{|y|, j/ E E \ a~ x E}, 

we have 7 > because (E \ a~ l E) C Ui 2 ^i(^), and £ Vi(-E). 2 < 2 < iV by [1, Lemma 5.23]. 
Recalling that rf/ = log iV/ log a we infer 

j<(\y\h) df , yeE„ j>2. (4) 

Consider rjj(z) = j E h(z,y) dfi(y), j > 1. For each n > we can take the representation 

Ej= U ^=^0-0^). 



l<fc<7V™ 

For 

Kk<N n 



the dominated convergence theorem [5j Proposition 7.1.1] implies that S^\z) A T)j(z), n — > 00, 
for each z G Z. 

For convenience, assume that diam.E = 1, then diamE^ = a~ n . If < < (5(h) — df/2, 
using ([3]) and the Cauchy-Schwarz inequality, we get 



\r]j(z)\ < \S m (z)\ +J2\ s{{n+m ( z ) -S (nj) (z)\ < \S (0 >\z)\ 

n>0 

n>0 l<fc<A r ™+ 1 

+**E E «^ (fc) lA*(^ iy) )l<l^.^)IIM^)l 

n>0 l<fc<A r ™+ 1 

+K,(^iV" +1 a 2 ^-^))) 1/2 (^a- 2 ^ £ |/i(4 (n+1)i) )| 2 ) V2 . (5) 

n>0 n>0 Kfc<iV n + 1 
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(The number k' is chosen such that E^ n+l ^' C E^'.) Lemma [T] implies 



a 

n>0 Kk<N n + 1 



oo a. s. 



Note that 2( / 9(/i) - /3) > d f = log iV/ log a. As in (jHJ) we obtain that 



n>m 



n>m 



n>m 



Kk<N n + 1 



consequently, the series 



n>0 



converges uniformly in z a.s. Each S^ n ^(z) has continuous paths, therefore r)j(z) has a continuous 
version. 



Recall that r)(z) = X^ ? >i r 7i( z )- Using 



and the Cauchy-Schwarz inequality we get for p > 1/2, 



h (z, 



< M for some constant M, estimate 



1/2 



IX*)i< £r 2 " 

j>m \j>m 
1/2 



Ef P M 2 

\j>m 



(EE E /'^K^ 

\j>m / \j>m n>0 l<fc<iV n + 1 



1/2 



1/2 



^((n+l)i) 



Now we apply Lemma [T] to the function sequences 

{gi(y), l>i} = {fMi E] ( y ), j>i}, 

{ 9l (y), 1>1} = {fa- nf3 l E (Xn +m (y), j > 1, n > 0, 1 < k < iV™ +1 } . 



(6) 



where Y^ii>i \gi(y)\ < C\y\ pdf (see (J4])). From the integrability assumption of our Lemma we infer 
that the all series in converge for m = 1. Thus we obtain that Y^j> m IVji^l tends to uniformly 
a.s. Therefore, i](z) has a continuous version. □ 



4 Stochastic heat equation on a nested fractal 

We consider the heat equation on E in the following mild sense 

u(t,x) = / p(t,x,y)u (y)dm(y) + ds p(t-s,x,y)f(s,y,u(s,y))dm(y) 

J E JO J E 

+ / d/i(y) / pit - s,x,y)a(s,y)ds, t G [0, T), xeE, (7) 

J E JO 

where u : [0, T] x E x Q — y R is an unknown measurable stochastic function, p is the transition 
density of the Brownian motion on E constructed in (4j, fi is a SM on 13(E) and uq is a random 
initial function. 
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We make the following assumptions. 

Assumption 2. uq : E x Q — > R is measurable and for each u G f2 the function mq(-,o;) is 
continuous and bounded on E, \uo(y, u)\ < C Uo (cu), and vanishes at oo. 

Assumption 3. / : [0, T] x E x R — > R is measurable and bounded, \f(s,y,r)\ < Cf. 
Assumption 4. f(s,y,r) is uniformly Lipschitz in y G E, r e R, 

yi, ri) - /(s,y 2 ,r 2 )| < - y 2 \ + \n -r 2 |) . 

Assumption 5. a : [0, T] x i£ — >■ R is measurable and bounded, |er(s, y) < Co- 
Assumption 6. a(s,y) is uniformly Holder continuous 'my E E, 

\a(s,y 1 )-a(s,y 2 )\ < K a \y x - y 2 f {a \ (3(a) > d f /2. 

Assumption 7. The spectral dimension d s of E is less than 4/3. 
Assumption 8. = a. s. for each x E E. 

Lemma 2. Under the Assumptions 1, 5, 6 and 7 the function 

H z ,v)= p(t- s,x,y)a(s,y)ds, z — (t, x) e [0, T] x E, y G E, 
Jo 

satisfies the conditions of Theorem^ for Z = [0, T] x E. 

Proof. We check the Holder continuity of h, the other conditions obviously hold true. By C we 
will denote a positive constant whose values may change. Using Assumption 6 and Theorem [T]b) 

we get 

\h(z,yi) - h(z,y 2 )\ < / \p(t-s,x,y 1 )-p(t-s,x,y 2 )\\<r(s,y 1 )\ds 

Jo 

+ / P{t ~ s,x,y 2 )\a(s,y 1 ) - a(s,y 2 )\ds 



o 



<C a \p(t - s,x,yx) -p(t - s,x,y 2 )\ds + CK G \y x - y 2 
Jo 

By Theorem [1] a) , for the first term we have 

nt-5 

\p(t-s,x,yx)-p(t-s,x,y 2 )\ds = 

Jo 

<C f \t- s)- 1 \y 1 -y 2 \ dw ~ df ds + C f [t - s)^' 2 ds 
Jo Jt-S 

< C(\ogT - log 5) \ yi - y 2 \ dw ~ df + C5 l - d °' 2 . 



Taking 8 = \yi — y 2 y w df ^^ 1 ds ^ 2 ' > and using that for all 7 > 0, t 7 logt — > as t — > 0, we obtain 
for any /3(h) < min{cL - d f , /3(a)}. If d s = 2d f /d w < 4/3 we can find (3(h) > df/2. □ 

Theorem 3. Suppose Assumptions 1-8 hold. 

1) Equation ffl) has a solution u(t,x). If v(t,x) is another solution, then for each t and x, 
u(t, x) = v(t, x) a. s. 

2) If, additionally, the function \y\ T be integrable w.r.t. \x on E for some r > df/2, then the 
solution u(t,x) has a version continuous on [0, T] x E. 
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Proof. Take u^°\t,x) = and set 

u {n+1) (t,x) = / p(t,x,y)u (y)dm(y) + I ds I p(t-s,x,y)f(s,y,u {n \s,y))dm(y) 



+ / dn{y) l p(t- s,x,y)a(s,y)ds, n>0. (8) 

JE Jo 

The measurability of the first and second summands of (jBJ) follows from the Fubini theorem. The 
third summand may be rewritten as limit in probability of the integrals of the simple functions 

/ hi(t,x,y)dfj,(y) , i ->■ oo, Kit, x, y) = ^ h k\t, a; ) 1 j4 w (v)- 
^ E k=i 



By (HI Lemma 2], the limit in probability of measurable processes is again measurable. 

Further, for all n, t, x we take the same version of the stochastic integral and obtain the 
following estimates for all u G Vt. Assumption 5 implies 

\u (n+1) (t,x) -u {n \t,x)\ 

= [ ds [_p{t-s,x,y)[f{s,y,u {n) {s,y))- /(s, u (n - 1} (s, y))] c/m(y) 

JO J E 

<K f ds p(t- s,x,y)\u in \s,y)-u ( - n ' 1 \s,y)\dm(y), n>2 
Jo Je 

Since / is bounded, we infer 

\u {2 \t,x) -u (1) {t,x)\ < 2C f [ ds fp(t-s,x,y)dm{y) = 2C f t. 

Jo Je 

Set g n (t) = sup x& j^\u^ n+1 \t, x) — u^ n \t,x)\, n>l. Then we get from (EJ), 



(9) 



9n(t)<K f g n ^i(s)ds, hence g n (t) < 2C f K n f - 



t 



n+l 



n + 1)! 



(10) 



i.e., the series Y^=o9 n ^J) converges uniformly on [0,T]. We define u(t,x) := lim^oo u^ n \t, x). 
Taking in (JSJ) the limit as n — > oo, we arrive at ([7]). 

If u(t,x) and v(t,x) both are solutions to ([7]), then we may repeat the arguments from Q 

to f fTUj) for u — v (instead of u( n+1 ) — u^) and g{t) = sup^.^ \u(t, x) — v(t, x)\ (instead of g n ) and 
obtain g = 0. 

From Lemma [2] and Theorem |2] it follows that has a continuous version. Then u(t, x) as 

uniform limit of u^ n \t, x) is continuous, too. □ 

Remark. For the Vicsek set Eys we have d s = log25/logl5 < 4/3, therefore Eys satisfies 
Assumption 7. It is known that Sierpifiski gasket is a nested fractal with d s = log9/log5 > 4/3, 
and Assumption 7 fails for this set. 
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